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SUMMARY 


A  computer  program,  SIMMIN,  which  Is  capable  of  finding  a  minimum  (If  one 
exists)  for  any  function  which  can  be  described  using  Fortran,  has  been  devel¬ 
oped  and  tested  on  an  IBM  PC  with  an  8087  math  coprocessor.  The  program  is  a 
modification  of  a  Pascal  algorithm  written  by  Caceci  and  Cacheris  and  modified 
by  Holloman  and  Otto.  A  simplex  method  determines  values  of  Initially  unknown 
coefficients  which  result  in  a  local  function  minimum. 

The  user  can  vary  initial  estimates  of  the  coefficients,  step  sizes,  max¬ 
imum  number  of  iterations,  and  the  convergence  criteria,  of  which  two  are 
availlable.  Final  values  of  the  coefficients,  an  estimate  of  the  function 
minimum,  and  the  number  of  iterations  required  are  displayed  on  screen,  while 
a  diagnostic  iteration  history  is  recorded  on  a  separate  file  which  can  be 
examined  after  each  execution. 

The  program  can  easily  be  converted  for  use  on  larger  computers,  perhaps 
allowing  the  use  of  single  precision.  Although  divergence  is  impossible,  a 
price  is  often  paid  in  more  iterations  and  computing  time  than  with  some  other 
traditional  minimization  methods. 
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I .  INTRODUCTION 

Techniques  for  minimizing  functions  include  stepwise  descent,  steepest 
descent,  Newton-Raphson  algorithms  (and  variations),  and  Marquardt  algorithms. 
However,  there  are  many  functions  for  which  these  well-known  methods  are  dif¬ 
ficult  to  implement,  are  computer  intensive,  and  are  prone  to  diverge.  Such 
functions  may  be  nonlinear,  partially  empirically  defined,  not  continuous  in 
regions  of  interest,  or  quite  complicated.  Inaccurate  calculation  of  deriva¬ 
tives  and  the  lack  of  good  initial  estimates  for  the  unknown  coefficients  are 
often  responsible  for  the  failure  of  these  methods. 

The  program  SIMMIN,  the  subject  of  this  report,  can  find  a  minimum  (if  one 
exists)  for  any  function  that  can  be  described  by  the  user  in  Fortran.  SIMMIN 
does  this  by  using  a  simplex  technique  to  calculate  values  of  initially 
unknown  coefficients  which  result  in  a  functional  minimum.  The  user  should 
have  some  familiarity  with  Fortran  so  that  the  subroutine  FUNC  in  the  source 
code  shown  in  the  Appendix  can  be  modified  for  the  function  of  interest.  A 
very  simple  function  has  been  included  for  instructional  purposes. 

II.  THE  MODEL 

The  simplex  methodology  for  finding  functional  minima  Is  used  in  SIMMIN. 
The  original  algorithm  was  written  by  Nelder  and  Mead  In  1965  [1]  and  was 
implemented  in  Extended  Mercury  Autocode  on  an  Orion  computer.  The  starting 
point  for  the  present  study  is  a  Pascal  algorithm  by  Caceci  and  Cacheris  to 
fit  curves  to  data  [2].  Holloman  and  Otto  modified  the  model  by  converting  to 
Fortran,  introducing  a  different  convergence  criterion,  to  be  discussed  later, 
and  hard-coding  user  inputs,  such  as  initial  estimates  and  step  sizes,  in 
order  to  make  the  model  more  transparent  to  the  user.  They  next  incorporated 
their  modified  model  Into  user-friendly,  curve-fitting  programs  for  IBM  PC 
[3,4].  In  addition  to  familiar  Fourier,  Lagrange  and  spline  methods,  data 
could  be  fit  to  any  equation  selected  by  the  user.  The  coefficients  in  the 
equation  were  determined  by  allowing  the  simplex  technique  to  minimize  the 
nonlinear  function  sum  of  squares  of  residuals  as  follows: 

SSR  =]C  t  Ycalc(X1,A1,A2,...,Amm)  ~  Yexper(Xj  )  ]2  (1) 

i  =  l 

where  Aj,  A2,...,Atnm  are  the  unknown  coefficients  and  N  is  the  number  of 
experimental  data  points,  Xj .  Least  squares  fitting  to  any  equation  was  done 
without  using  derivatives  or  matrix  operations.  However,  it  was  noted  that 
for  most  fitting  functions,  execution  times  Increased  greatly  with  the  number 
of  coefficients,  above  approximately  three. 

In  the  present  study,  the  fitting  aspects  of  the  model  were  removed.  User 
control  over  most  inputs  was  restored  since  finding  a  minimum  (or  minima) 
usually  requires  experimentation  with  Initial  estimates,  step  sizes,  the  maxi¬ 
mum  number  of  iterations,  and  the  convergence  criteria,  of  which  two  are 
available. 

The  user  is  referred  to  References  1  and  2  for  details  of  the  model,  but  a 
short  introduction  is  offered  here.  A  simplex  is  a  geometric  figure  which  has 
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one  more  vertex  than  the  space  in  which  it  is  defined  has  dimensions.  For 
example,  the  user  will  find  in  the  Appendix  (the  source  code  for  program 
SIMMIN)  a  function  of  the  following  form  defined  in  line  numbers  281-284: 

F  =  (A-l.)2  +  (B-2.)2  +  (C-3.)2  (2) 

The  simplex  for  this  function  would  have  four  vertices  in  three-dimensional 
space.  Each  vertex  is  characterized  by  four  values:  A,B,C,F.  Subroutine 
SIMPLEX  defines  the  simplex  as  a  two-dimensional  array  SIMP(I,J),  with  I  being 
the  vertex  index  and  J  the  dimension  index.  It  should  be  clear  that  for  our 
example,  values  A  =  1,  B=2,  C  =  3  yield  the  one  and  only  function  minimum  F 
=  0. 

The  principal  of  the  simplex  method  is  to  move  and  reshape  the  simplex  in 
ways  resulting  in  decreased  values  of  the  response  function  as  iterations 
proceed.  Each  iteration  involves  combinations  of  reflection,  expansion, 
contraction,  and  shrinkage  of  the  vertices  according  to  logical  rules,  and 
then  testing  the  newly  calculated  function  value  against  the  previous  value. 

If  the  new  value  is  better  (meaning  lower)  than  the  previous  one,  then  the 
vertices  are  accepted  and  the  simplex  proceeds  in  the  same  general  direction. 
If  not,  then  the  vertex  changing  logic  proceeds  in  an  orderly,  efficient 
search  or  a  new  direction. 

III.  COMPUTER  IMPLEMENTATION 

The  present  study  was  done  on  an  IBM  PC  equipped  with  an  8087  coprocessor 
using  Ryan-McFarland  Fortran  (V2.ll).  This  Fortran  creates  very  efficient 
executable  code  at  the  expense  of  longer  compilation  times.  Double  precision 
was  used  for  all  reals  by  appropriately  placed  implicit  double  precision 
statements.  The  user  can  easily  remove  the  double  precision  If  SIMMIN  is  used 
on  fast  computers  with  large  word  size,  where  large  numbers  of  Iterations  and 
truncation/round-off  errors  are  not  likely  to  cause  problems  using  single 
precision. 

The  primary  modification  the  user  will  make  is  to  substitute  a  function 
of  interest  for  the  demonstration  function  in  subroutine  FUNC  (see  the 
Appendix,  Lines  272-287).  The  program  requires  that  the  coefficients  enter 
the  subroutine  as  an  array  COEF(I),  I  *  1,...,MM,  where  MM  is  the  number  of 
coefficients  to  be  determined.  The  user  is  then  free  to  convert  to  a  more 
convenient  notation.  Perhaps  the  function  to  be  minimized  is  available  as  a 
function  subprogram  or  a  subroutine  which  would  then  be  called  from  within 
subroutine  FUNC.  The  calculated  value  of  the  function  must  be  returned  from 
FUNC  as  C0EF(MM+1). 

Constraints  can  be  added  in  the  form  of  Fortran  statements  directing  the 
program  flow  to  statement  999  in  FUNC,  often  In  the  form  of  IF  statements 
restricting  the  coefficients  or  function  to  certain  regions.  Note  what  hap¬ 
pens  then:  the  function  value  C0EF(MM+1)  is  returned  as  the  previous  value 
plus  a  very  large  number  1.D30.  Most  calls  of  subroutine  FUNC  from  within 
subroutine  SIMPLEX  are  followed  by  IF  tests  to  determine  if  the  new  trial 
function  value  RNEXT(MM+1)  is  less  than  the  previous  value,  referred  to  in 
subroutine  simplex  as  SIMP(MM+1 , IH(MM+1) ) .  It  certainly  will  not  be  if  a 
constraint  condition  has  been  exceeded  and  the  logic  will  guide  the  simplex 
away  from  this  region. 
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An  important  additional  word  about  constraints:  the  user  may  well  have 
them  without  being  aware  of  them  for  awhile.  The  simplex  moves  about  in 
hyperspace  without  regard  to  such  problems  as  attempting  to  take  logarithms  of 
zero  or  of  negative  numbers,  dividing  by  zero,  or  attempting  to  evaluate  func¬ 
tions  (such  as  factorials,  gamma  functions,  etc.)  in  regions  leading  to  com¬ 
puter  underflow  or  overflow.  The  user  should  protect  against  all  such  possi¬ 
bilities  by  adding  appropriate  statements  to  direct  the  program  to  statement 
999,  or  to  do  the  equivalent  action  described  after  statement  999  somewhere 
else. 

A.  Input  Parameters 

Values  supplied  by  the  user  during  interactive  operation  of  SIMMIN 
are  shown  within  boxes  in  Figures  1,  2,  and  3. 

1.  MM  -  The  number  of  coefficients  to  be  calculated,  3  for  the 
example  in  Equation  (2). 

2.  MAXITR  -  The  maximum  number  of  iterations  allowed.  This  may  be 
only  approximate  since  the  program  now  counts  every  trial  set  of  vertices  as 
an  iteration,  regardless  of  whether  the  function  value  increases  or  decreases. 
Therefore,  the  number  of  iterations  shown  as  part  of  the  results  could 
slightly  exceed  MAXITR. 

3.  Initial  Estimates  of  the  Coefficients  -  The  user  must  insure  that 
their  order  conforms  to  the  definition  of  coefficients  in  subroutine  FUNC. 

One  should  always  experiment  with  these  initial  estimates  to  see  the  effect  on 
results . 

4.  Convergence  Criterion  Option 

ICON  =1  —  The  program  compares  relative  errors  of  calculated 
coefficients  with  a  limit  entered  by  the  user,  typically  between  l.D-05  and 
l.D-04.  If  every  relative  error  is  less  than  the  limit,  then  convergence  is 
obtained.  Before  defining  this  error,  recall  for  Equation  (2)  that  there  are 
four  vertices  for  the  simplex:  A,  B,  C,  F,  and  that  each  vertex  is  charac¬ 
terized  by  four  values:  A,  B,  C,  F.  The  simplex  is  defined  as  a  two- 
dimensional  array  SIMP(I,J),  with  I  being  the  vertex  index  and  J  being  the 
dimension  index.  The  error  is  defined  as  follows  (Appendix,  Line  231): 

ERR0R( I)=(SIMP( I , IH( I) )-SIMP( I ,L( I) ) )/SIMP( I , IH( I) )  (3) 

subroutine  ORDER  has  provided  IH(I)  and  L(I),  the  indices  of  the  highest  and 
lowest  values,  respectively,  currently  associated  with  vertex  I.  The  value  of 
ERROR(MMfl),  the  function  error,  does  not  seem  to  decrease  to  nearly  as  small 
as  the  coefficient  errors  ERROR(I),  1=1, MM  even  though  Inspection  of  latter 
iterations  in  file  ITER. OUT  shows  a  minimum  is  close  to  being  found,  by  either 
convergence  criterion.  Therefore,  the  error  comparisons  were  restricted  to 
the  coefficients  only. 

ICON  =2  —  This  convergence  option  allows  the  program  to  keep  a 
count  of  the  number  of  consecutive  times  that  new  trial  vertices  have  been 
calculated  without  reducing  the  value  of  the  function.  If  this  count  exceeds 
(MM+1)*3  then  convergence  is  obtained.  It  was  arrived  at  through  experience 
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with  nonlinear  functions  in  fitting  applications,  and  relieves  the  user  from 
choosing  acceptable  limits  for  errors  [3,4].  The  simplex  has  usually  found  a 
minimum.  However,  it  may  be  trapped  near  some  unusual  topological  feature, 
perhaps  moving  up  and  down  a  valley  where  coefficient  values  change  but  the 
function  value  does  not.  When  dealing  with  new,  unfamiliar,  complicated 
functions,  it  is  always  advisable  to  experiment  with  initial  estimates  and 
step  sizes  and  to  examine  at  least  the  latter  iterations  in  file  ITER. OUT. 

B.  Results 

In  Figures  1,  2,  and  3,  test  results  begin  after  the  lines 
"*****CALCULATING  SIMPLEX*****".  Tables  1  and  2  show  the  latter  iterations 
from  iteration  history  file  ITER. OUT  produced  by  executing  the  cases  shown  in 
Figures  1  and  2,  respectively.  After  each  execution,  this  file  can  be 
examined,  printed,  or  renamed,  but  file  ITER. OUT  is  destroyed  and  recreated 
near  the  beginning  of  each  execution. 

Figure  1  shows  convergence  criterion  1  being  satisfied  with  74  itera¬ 
tions  for  a  selected  error  limit  of  l.D-05.  Table  1  shows  that  all  coef¬ 
ficient  relative  errors  are  indeed  less  than  l.D-05  at  iteration  74.  This 
table  also. shows  the  construction  of  this  diagnostic  file,  with  results 
labeled  for  iteration  number  69. 

Figure  2  shows  the  same  case  as  Figure  1  except  that  convergence  cri¬ 
terion  2  was  selected.  Note  that  the  final  value  of  the  function 
(0.414162E-24)  is  lower,  and  the  number  of  iterations  (224)  is  higher. 
Criterion  2  is  often  more  stringent  than  criterion  1. 

Table  2  shows  iterations  211-224  from  file  ITER. OUT.  Notice  that 
beginning  with  iteration  213,  the  relative  errors  develop  a  pattern  as  itera¬ 
tions  proceed:  the  first  and  third  ones  switch  back  and  forth  between  pairs 
of  values,  and  the  second  remains  constant.  This  continues  for  3*(MM+1)=*12 
iterations,  and  then  criterion  2  Is  satisfied.  This  type  of  error  behavior  is 
common  when  the  simplex  has  locked  in  to  a  minimum.  Given  the  logical  rules 
for  changing  vertex  values  and  the  truncation/round-off  limits  in  the  com¬ 
puter's  arithmetics,  the  program  is  unable  to  find  a  change  which  results  in  a 
lower  function  value.  Without  the  logic  for  criterion  2,  iterations  would 
have  proceeded  until  the  iteration  limit  was  exceeded. 

The  case  in  Figure  3  is  the  same  as  that  in  Figure  2  except  the  Ini¬ 
tial  estimates  were  200,  300,  400  instead  of  2,  3,  4,  which  were  much  closer 
to  the  correct  answers  1,  2,  3.  Using  convergence  criterion  2,  268  iterations 
were  required.  Approximate  wall  clock  execution  time  was  30  seconds  for 
Figure  3  and  26  seconds  for  Figure  2  (224  iterations).  For  comparison,  using 
single  precision,  the  execution  times  were  approximately  16  and  12  seconds, 
respectively.  However,  the  estimation  of  the  true  function  minimum  of  0.0 
was  about  10-^2  -  10^  using  single  precision,  while  double  precision  improved 
it  to  approximately  102^  -  10^0. 

IV.  CONCLUSIONS 

The  computer  program  SIMMIN  is  capable  of  finding  a  minimum  (If  one 
exists)  for  any  function  which  can  be  described  using  Fortran.  It  was  tested 
on  an  IBM  PC  with  an  8087  math  coprocessor  using  Ryan-McFarland  Fortran.  This 
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simplex  method  determines  values  of  initially  unknown  coefficients  which 
result  in  a  local  function  minimum. 

The  user  can  vary  initial  estimates  of  the  coefficients,  step  sizes, 
maximum  number  of  iterations,  and  the  convergence  criteria,  of  which  two  are 
available.  Final  values  for  the  coefficients,  an  estimate  of  the  function 
minimum,  and  the  number  of  iterations  required  are  displayed  on  the  screen.  A 
diagnostic  iteration  history  is  recorded  on  a  separate  file  which  can  be  exa¬ 
mined  after  each  execution. 

The  program  can  easily  be  converted  for  use  on  larger  computers  with 
larger  word  sizes,  perhaps  allowing  the  use  of  single  precision.  Although 
divergence  is  impossible  using  this  simplex  method,  a  price  is  often  paid  in 
more  iterations  and  computing  time  than  with  some  other  traditional  minimiza¬ 
tion  methods,  such  as  those  mentioned  in  the  introduction. 

The  user  is  cautioned  against  drawing  overly  optimistic  conclusions  from 
the  test  results.  A  very  simple  test  function  was  used.  Actual  functions 
encountered  in  research  and  development  -  analytical,  empirical,  or  hybrids  - 
which  cause  problems  when  using  minimization  methods  such  as  those  mentioned 
in  the  introduction  could  require  care  when  using  simplex  techniques  as  well. 

Here  are  a  few  advantages  simplex  methods  such  as  those  in  SIMMIN  offer 
over  other  methods : 

1.  No  derivatives  of  any  kind  are  calculated. 

2.  Divergence  is  impossible. 

3.  Initial  estimates  for  the  coefficients  can  be  poor,  but  this  will 
slow  down  convergence,  depending  upon  step  size. 

4.  Simplex  logic  insures  economical  calculation  of  the  response 
function. 

5.  No  matrix  operations  are  involved. 

6.  Constraint  conditions  which  can  be  stated  in  Fortran  can  be  added 
very  easily. 

Here  are  a  few  disadvantages  of  simplex  compared  with  other  methods: 

1.  More  iterations  and  computing  time  can  be  required.  SIMMIN  is 
recommended  for  larger  computers,  especially  when  the  function  is 
complicated,  when  other  methods  have  proved  unsatisfactory  or  need 
to  be  verified,  or  when  the  user  has  little  knowledge  of  the  func¬ 
tional  behavior. 

2.  The  following  can  cause  computational  problems:  very  large  initial 
estimates  coupled  with  small  steps,  insensitivity  of  the  function 
to  one  or  more  coefficients,  too  large  an  acceptable  error  when 
using  an  error  comparison  convergence  criterion.  These  and  other 
problem  areas  are  discussed  in  much  more  detail  in  Reference  2. 
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ENTER  MM,  NUMBER  OF  PARAMETERS  TO  BE  CALCULATED  -131 

ENTER  MAXITR,  MAXIMUM  NO.  OF  ITERATIONS,  TYPICALLY  SEVERAL  HUNDRED.  PROGRAM 
USUALLY  CONVERGES  WITHIN  20*(MM**2)  ITERATIONS,  -flffffl 

ENTER  MM  INITIAL  ESTIMATES  OF  PARAMETERS.  BE  SURE  THE  ENTRY  ORDER  CONFORMS 
WITH  YOUR  DEFINITION  OF  PARAMETERS  COEF(I)  IN  SUBROUTINE  FUNC.  ESTIMATES 
USUALLY  CAN  BE  QUITE  INACCURATE.  - 

[273741 


ENTER  MULTIPLIER  IN  STEP  CALCS.,  APPROX.  RANGE  FROM  0.1  to  0.5. 
ESTIMATE  *  MULTIPLIER)  -  [7T] 


(STEP 


ENTER  CONVERGENCE  CRITERION  OPTION  - 

1,  FOR  COMPARISON  OF  RELATIVE  ERRORS  WITH  INPUT  LIMIT. 

2,  FOR  THE  CONDITION  THAT  THE  NUMBER  OF  UNSUCCESSFUL  ATTEMPTS  TO  DECREASE  THE 
FUNCTION  EXCEEDS  3*(MM+1)  -(T| 

ENTER  LIMIT  ON  RELATIVE  ERRORS,  THE  SAME  FOR  ALL  COEFFICIENTS. 

RECOMMENDED  RANGE  l.D-06  to  l.D-04.  -  |l.D-0~5) 

*****  CALCULATING  SIMPLEX***** 

FINAL  VALUES  OF  COEFFICIENTS  - 
1.00000  2.00000  3.00000 

FINAL  VALUE  (MINIMUM)  OF  FUNCTION  =  0.417040E-10 
NUMBER  OF  ITERATIONS  =74 
Execution  terminated  :  OK 


1 


Figure  1.  Screen  input  and  results  for  a  typical  execution  of  SIMMIN. 
Convergence  criterion  1  was  selected. 


w 


>  w*.*  .  m\  •>  v,"  »' 


*  V  Vfc* 


TABLE  1.  Latter  Contents  of  Iteration  History  File  ITER. OUT  Produced  During 
Execution  With  Input  Parameters  Shown  in  Figure  1. 

69  =  Iteration  number 
1.000 

2.000  Current  values  of  coefficients 
3.000 

0  l?10E-09  Estimate  of  function  value,  printed  only  when  the  iteration  led 
to  a  better  (lower)  function  value. 

1 .48453975763265370E-05  7.86123165807031016E-06  7.56428230863572653E-06 

70  \  / 

Estimates  of  relative  errors  for  the  coefficients,  MM  =  3 

1.0000 

2.000 

3.000 

0.7947E-10 

1 . 48453976763265370E-05  4 . 89503933556161 228E-06  7 . 56428230863572653E-06 

71 

1.000 

2.000 

3.000 

1.44438985795023159E-05  4.89503933556161228E-06  7 .56428230863572653E-06 

72 

1.0000 

2.000 

3.000 

0.3374E-10 

1 .20481 762607656796E-05  3. 1 281 55729957 17666E-06  5.20121882225269973E-06 

73 

1.000 

2.000 

3.000 

0. 3160E-10 

1 .20481762607656796E-05  3 . 1 281 55729957 17666E-06  2 . 91527801458411470E-06 

74 


1.000 

2.000 

3.000 


All  relative  errors  <  l.D-05. 

Option  1  convergence  (ICON  =  1)  obtained. 


6.9  3297636726561 162F.-06  3 . 1 28 1 55729957 1 7666E-06  2 . 9152780145841 1470E-06 


- >2 
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ENTER  MM,  NUMBER  OF  PARAMETERS  TO  BE  CALCULATED  -[3J 

ENTER  MAXITR,  MAXIMUM  NO.  OF  ITERATIONS,  TYPICALLY  SEVERAL  HUNDRED.  PROGRAM 
USUALLY  CONVERGES  WITHIN  20*(MM**2)  ITERATIONS.  -  |300{ 


ENTER  MM  INITIAL  ESTIMATES  OF  PARAMETERS.  BE  SURE  THE  ENTRY  ORDER  CONFORMS 
WITH  YOUR  DEFINITION  OF  PARAMETERS  COEF(I)  IN  SUBROUTINE  FUNC.  ESTIMATES 
USUALLY  CAN  BE  QUITE  INACCURATE.  - 


APPROX.  RANGE  FROM  O.I  to  0.5.  (STEP 


ENTER  MULTIPLIER  IN  STEP  CALCS. 
ESTIMATE  *  MULTIPLIER)  -  [72 


ENTER  CONVERGENCE  CRITERION  OPTION  - 

1,  FOR  COMPARISON  OF  RELATIVE  ERRORS  WITH  INPUT  LIMIT. 

2,  FOR  THE  CONDITION  THAT  THE  NUMBER  OF  UNSUCCESSFUL  ATTEMPTS  TO  DECREASE  THE 
FUNCTION  EXCEEDS  3*(MM+1)  -  jT] 

*****  CALCULATING  SIMPLEX***** 

FINAL  VALUES  OF  COEFFICIENTS  - 
1.00000  2.00000  3.00000 

FINAL  VALUE  (MINIMUM)  OF  FUNCTION  =  0.414162E-29 
NUMBER  OF  ITERATIONS  =  224 
Execution  terminated  :  OK 


Figure  2.  Screen  input  and  results  for  a  typical  execution  of  SIMMIN. 
Convergence  criterion  2  was  selected. 


yy.yv/y-'-y' 
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a 
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TABLE  2.  Latter  Contents  of  File  ITER. OUT  Produced  by  Input  Parameters  in 
Figure  2.  Repetitive  Pattern  in  Coefficient  Errors  Begins  at 
Iteration  213. 


1.000 

2.000 

3.000 

0.4930E-31 


2. 22044604923031 111E-1 5 
212 


1.11022302462515536E-15 


8 . 881 784 197001 24936E-16 


1.0000 

2.000 

3.000 


2. 22044604925031 111E-1 5 
213 


6.66133814775093530E-16 


5. 92 1 189464667 501 2 2E- 16 


1.0000 

2.000 

3.000 


1.99840144432528138E-15 

214 


6 . 66133814775093530E-16 


5. 92 11 89464667 5002 3E-1 6 


1.000 

2.000 

3.000 


2 . 22044604925031 11 IE-15 
215 


6.66133814776093530E-16 


5. 921 189464667501 22E-16 


1.0000 

2.000 

3.000 


1 . 99840 144432 528 138E- 15 
216 


6.661338147750935  30E- 1 6 


5.92118946466750023E-16 


1.000 

2.000 

3.000 


2. 22044604825031 111E-15 
217 


6.66133814775093530E-16 


5.92118946466750122E-16 


1.0000 

2.000 

3.000 


I 


I 


1 • 99840144432528138E-15 
218 

1.000 

2.000 

3.000 

2 .22044604925031 11 IE-1 5 

219 

1.0000 

2.000 

3.000 

1.99840144432528138E-15 

220 

1.000 

2.000 

3.000 

2. 22044604925031 111E-15 
221 

1.0000 

2.000 

3.000 

1.99840144432528138E-15 

222 


TABLE  2.  (Concluded) 
6.66133814775093530E-16 


6.66133814775093530E-16 


6.66133814775093530E-16 


6.66133814775093530E-16 


6. 6613381477 5093530E-16 


5. 921 189464667 50023E-1 6 


5. 921 189464667 501 22E-1 6 


5.92118946466750023E-16 


5. 921 189464667 501 22E-1 6 


5.92118946466750023E-16 


1.000 

2.000 

3.000 


2. 22044604925031 111E-15  6.66133814775093530E-16  5.92118946466750122E-16 

223 


1.0000 

2.000 

3.000 

1.99840144432528138E-15  6.66133814775093530E-16  5.92118946466750023E-16 

224 

1.000 

2.000 

3.000 


2. 22044604925031 111E-1 5 


6.66133814775093530E-16 


5.92118946466750122E-16 


ENTER  MM,  NUMBER  OF  PARAMETERS  TO  BE  CALCULATED  ~(T| 

ENTER  MAXITR,  MAXIMUM  NO.  OF  ITERATIONS,  TYPICALLY  SEVERAL  HUNDRED. 
PROGRAM  USUALLY  CONVERGES  WITHIN  20*(MM**2)  ITERATIONS.  -  (5551 


ENTER  MM  INITIAL  ESTIMATES  OF  PARAMETERS.  BE  SURE  THE  ENTRY  ORDER  CONFORMS 
WITH  YOUR  DEFINITION  OF  PARAMETERS  COEF(I)  IN  SUBROUTINE  FUNC.  ESTIMATES 
USUALLY  CAN  BE  QUITE  INACCURATE.  - 
[200, 30Q,400| 

ENTER  MULTIPLIER  IN  STEP  CALCS.,  APPROX.  RANGE  FROM  0.1  to  0.5.  (STEP  - 
ESTIMATE  *  MULTIPLIER)  -  [7T| 

ENTER  CONVERGENCE  CRITERION  OPTION  - 

1,  FOR  COMPARISON  OF  RELATIVE  ERRORS  WITH  INPUT  LIMIT. 

2,  FOR  THE  CONDITION  THAT  THE  NUMBER  OF  UNSUCCESSFUL  ATTEMPTS  TO  DECREASE  THE 
FUNCTION  EXCEEDS  3*(MM+1)  -JT] 

*****  CALCULATING  SIMPLEX***** 

FINAL  VALUES  OF  COEFFICIENTS  - 
1.00000  2.00000  3.00000 

FINAL  VALUE  (MINIMUM)  OF  FUNCTION  =  0.246519E-30 
NUMBER  OF  ITERATIONS  =  268 
Execution  terminated  :  OK 


Figure  3. 


Screen  input  and  results  for  a  typical  execution  of  SIMMIN. 
Initial  coefficient  estimates  were  100,  200,  300. 
Convergence  criterion  2  was  selected. 


APPENDIX 

FORTRAN  SOURCE  CODE  FOR  SIMMIN 


) 

I 


a 


A 


s 


A-l/(A-2  Blank) 


\ 

! 


r  V> 


iij 


1 


PROGRAM  SIMMIN 

FROGRAM  3 1  MM  IN  CAN  FIND  A  rffNIMCM  <  IF' ONE  EXISTS)  FOR  ANY 
FUNCTION  DEFINED  ST  THE  USER  IN  SUBROUTINE  FUNC .  THE  FUNCTION 
CAN  BE  LINEAR  OR  NONLINEAR,  PARTIALLY  EMPIRICALLY  DEFINED,  AND 
NEED  NOT  3E  CONTINUOUS.  SIMMIN  OOES  THIS  BY  USING  SIMPLEX 
TECHNIQUES  TO  CALCULATE  THE  UALUES  OF  INITIALLY  UNKNOWN 
COEFFICIENTS  THAT  YIELD  THE  MINIMUM.  PROGRAM  IS  BASED  ON  A 
PASCAL  ALGORITHM  IN  "FITTING  CURVES  TO  DATA",  BYTE,  MAY  1664, 

WHICH  WAS  CONVERTED  Tu  FORTRAN??  AND  MODIFIED  FOR  USE  IN 
GENERAL  PC  CURVE  FITTING  PROGRAMS  BY  MILES  HOLLOMAN  AND  BILL  OTTO, 
DIRECTED  ENERGY  DIRECTORATE,  MICOM.  FURTHER  MODIFICATIONS 
RESULTED  IN  THIS  INITIAL  VERSION  OF  SIMMIN  BY  JERRY  D.  SMITH, 

OF  THE  SAME  DIRECTORATE.  USER  SHOULD  BE  FAMILIAR  ENOUGH  WITH 
FORTRAN  TO  MODIFY  SUBROUTINE  FUNC  AS  REQUIRED ,  CONVERGENCE 
CRITERION-  IF  DESIRED  (SEE  COMMENTS  WITHIN  SUBROUTINE  SIMPLX), 

AND  PERHAFS  MINOR  FORMAT  CHANGES  IF  MANY  UNKNOWN  COEFFICIENTS 
COEF(I)  OCCUR  IN  THE  FUNCTION.  CURRENT  MAX.  NO.  =  10. 

NOTE  -  EXECUTION  TIMES  INCREASE  GREATLY  WITH  NUMBER  OF 

COEFFICIENTS  ,  ABOVE  ThuEE. 

AFTER  EXECUTION  FILE  ITER. OUT  WILL  CONTAIN  A  HISTORY  OF  THE 
ITERATION  RESULTS  :  COEFFICIENT  AND  FUNCTION  VALUES,  AND  RELATIVE 
ERROR  VALUES  FOR  THE  COEFFICIENTS.  USER  CAN  EXAMINE  AFTER  EXECU¬ 
TION  IF  DESIRED.  IT  IS  DESTROYED  AND  RECREATED  DURING  THE  NEXT 
EXECUTION. 

ADVANTAGES  OVER  OTHER  NONLINEAR  METHODS: 

1.  NO  DERIVATIVES  ARE  CALCULATED. 

2.  DIVERGENCE  IS  IMPOSSIBLE. 

3.  INITIAL  ESTIMATES  FOR  COEFFICIENTS  CAN  BE  QUITE  POOR, 

ALTHOUGH  THIS  WILL  SLOW  OOUN  CONVERGENCE. 

4.  SIMPLEX  LOGIC  INSURES  ECONOMICAL  CALCULATIONS  OF  RESPONSE 
FUNCTION. 

5.  NO  MATRIX  OPERATION  IS  INVOLVED. 

5.  CONSTRAINTS  CAN  BE  ADDED  VERY  EASILY.  ( SEE  SUBROUTINE  FUNC.' 
DISADVANTAGES  OVER  OTHER  NONLINEAR  METHODS: 

!  .  MORE  ITERATIONS  AMD  COMPUTING  TIME  ARE  OFTEN  REQUIRED.  SIMMIN 
RECOMMENDED  FOR  FAST  PC’S,  AMD  NUMBER-CRUNCHERS ,  ESPECIALLY 
WHEN  THE  FUNCTION  TO  EE  MINIMIZED  IS  COMPLICATED  AND  USER  HAS 
VERY  LITTLE  IF  ANY  KNOWLEDGE  OF  THE  FUNCTIONAL  BEHAVIOR. 

2.  FOLLOWING  CAN  CAUSE  COMPUTATIONAL  PROBLEMS,  AND  RESULT  IN 
INACCURATE  ANSWERS:  VERY  LARGE  INITIAL  GUESSES,  VERY  SMALL 
INITIAL  INCREMENTS,  INSENSITIVITY  OF  FUNCTION  TO  ONE  uR  MORE 
COEFFICIENT* ,  TOG  LARGE  TN  ACCEPTABLE  ERROR  i WHEN  USING 
CONVERGENCE  OPTION  ICON-2). 


JERRY  0.  SMITH 
! I  MAY  e7 


COMMON  FINAL!  !  1  ) ,SIMP( 11,11)  ,STEP(  1  I  ) ,MM ,LUO .MAXITR  .NITER  .ICON  ,RX 
IMPLICIT  DOUBLE  PRECISION  (A-H.0-2) 


P 


■'-.v --- --V  ZK'C'T 


A-3 


yrz^v  Tijmrx  *xrxr  nr  wj- 


rowuflmviw.  r  w  *? 


> » n.'i w.k1  w  wuvmu.' w  P-V-vm  ■.“  y  j.vj  ^v-vjvjw  vj  yj  ww 


57 

58 

59 

60 

SI  C 
62 

S3  I SB3 

64 

65 

SB  1601 
S7 
68 

69 

70 

71 

72  1602 

73 

74 

75 

76 

77 

73  1603 
75 
80 
81 

82  13 

83  C 

84 

85  1610 

86 

87 

88 
83 

90 

91 

32  1612 

33 

34 

95 

96  C 

97 

53  C 
53 

100  1504 

101 

1C2  1606 

!  0  j 

1  04  1 b  0  d 

105  C 

106 
107 


LUO-7 

OPEN (UNIT -LUO, FILE-’ ITER. OUT’ .STATUS- ' UNKNOWN '  > 

CLOSE! UNIT-LUD. STATUS- ’DELETE’ ) 

0PEN( UNIT-LUQ  .PILE-’ ITER. OUT' .STATUS-'NEW’ ) 

WRITE!*  ,1600) 

FORMAT* IHt  ,  ’ENTER  MM,  NUMBER  OF  PARAMETERS  TO  BE  CALCULATED  -  ') 
READ! •  ,«  )  MM 
WRITE! *  ,1601 ) 

FORMAT! IX  ,  'ENTER  MAXITfi,  MAXIMUM  NO,  OF  ITERATIONS  .TYPICALLY  SEVER 
XAL  HUNOREO.’,/’  PROGRAM  USUALLY  CONVERGES  WITHIN  20*!MMW)  ITERAT 
XICNS.  -  '  ) 

READ**,*!  M.TXITR 
MM1-MM+1 
WRITE! »,  1 602) 

FORMAT! IX  , ’ENTER  MM  INITIAL  ESTIMATES  OF  PARAMETERS.  BE  SURE’ 

X,/,'  THE  ENTRY  ORDER  CONFORMS  WITH  YOUR  DEFINITION  OF  PARAMETERS' 
COEFili  IN  SUBROUTINE  FUND.  ESTIMATES  USUALLY  CAN  BE  QUITE' 
INACCURATE  -  ’/> 

READ!  *  ,•  )  (SIMP!  I  ,1  )  ,1-1  ,  MM  > 

WRITE. ♦ ,1603) 

FORMAT (IX, ’ ENTER  MULTIPLIER  IN  STEP  CALCS.,  APPROX. RANGE  FROM  0.1 
a  ’  /,’  TO  0-5.  (STEP-ESTIMATE  *  MULTIPLIER )  -  ’> 

READ! »  ,* )  SMULT 

DO  10  J  =  1  .Mil 

3TEP<  J  )-SIMF( J  ,1  )* SMULT 


PRINT  1610 

FORMAT! IX,  'ENTER  CONVERGENCE  CRITERION  OPTION  -’,/1X,  M,  FOR  COM 
XpARISON  OF  RELATIVE  ERRORS  WITH  INPUT  LIMIT.’, /IX,  *2,  FOR  THE  CON 
XOITION  that  the  NUMBER  OF  UNSUCCESSFUL  ATTEMPTS  TO', /IX,  ’  DECRE 
XASE  THE  FUNCTION  EXCEEDS  3*<MM+!  >  -  ’> 

READ! *  ,*  >  ICON 
IF! ICOK.EQ. 1 i  THEN 
PRINT  1612 

FORMAT! IX,  'ENTER  LIMIT  ON  RELATIVE  ERRORS,  THE  SAME  FOR  ALL  COEFF 
XICIENTS.  ’  ./IK.  'RECOMMENDED  RANGE  5 .D-0S  TO  1 .D-04 
READ! * , *  >  RX 
END  IF 

CALL  SIMPLEX 

WRITE'.  »  ,  1604  )  (FINAL!  I  ?  ,1-1  ,M11 ) 

FORMAT! IX  , 'FINAL  VALUES  OF  COEFFICIENTS  -  '  ,/3X  ,  6612.65 
PRINT  1 S06  ,  FINAL! MM1  > 

FORMAT!  '  FINAL  VALUE  (MINI MUM  >  OF  FUNCTION  -  ’,  £12.5.' 

PRINT  t  S08  ,  til  IER 

FORMAT! IX  ,  ’NUMBER  OF  ITERATIONS-’ . U > 


STOP  ’OK’ 
END 


4 


5 

I 


% 


SUBROUTINE  SIMPLEX 


108 
102 
110 
1  1  1 
112 

113 

114 
1  IS 
I  tS 
1  17 
119 

119 

120 
121 
122 

123 

124 

125 
I2S 

127 

128 
129 
1 3a 
131 
t32 
1 33 


1030 

C 


1040 

C 

C 

c 

c 


LOGICAL  DONE .LEO 

COMMON  FINAL!  1  1  )  ,SIMP<  11,1!)  .STEP!  I  1  >  ,1411  .LUO  .I1AXITR  ..NITER  ,  ICON  ,RX 
DIMENSION  CENTER! 1 1  •  .ERROR! 1 1  i  .RMAXERi I  1  )  ,PS< 1 1)  ,Q< 1 1 > .  IH( 1  I  ) 
DIMENSION  RNEXT! 1 1 ).L( 11 ) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

OLOU-1 .30+38 
OLQWL’I .00+38 
LCNT-0 
PRINT  1030 

FORMAT!’  *•»»*  CALCULATING  SIMPLEX  ’> 

MM1-MM+! 

CO  1040  1  =  1  ,1  I 
RMAXER! I >=RX 

FOLLOWING  ARE  PARAMETERS  TO  UARY  THE  SIMPLEX  MODIFICATIONS  AND 
SHOULD  GENERALLY  SE  ACCEPTED  AS  OPTIMAL  UALUES. 

ALFA  -1.0 
SETA  =  0.5 
GAMMA  =2.3 


1 34 

TUO  =  2.0 

135 

root:  -  DSQRT! TWO  > 

133 

CALL  FUND!  SIMP'  1.1)  .MM1  ) 

137 

00  1062  1=1,  MM 

138 

RNN  »  MM  1 

139 

PS!  I  )*STEP-'  I  >  +  <  <S<3RT<FNN  1+MM-! 

140 

0!  I  )=STEP(  I  )•<  (SORT!  RNN)  >-l  >/< 

141 

I860 

CONTINUE 

142 

DO  1063  1=2 ,MMI 

143 

DO  1070  1= 1  ,MM 

1  44 

1070 

SIMP! J ,1  )=SIMP' J  .  1  >+0<  J  ) 

145 

SIMP!  1-1,1  >=SIMP<  1-1,1  )+PS!  1-1 

US 

CALL  FUNC! SIMP! 1,1) ,MMI ) 

i  47 

1033 

CONTINUE 

!  43 

OC  1090  1=1 ,MMt 

149 

LC I  )  =  1 

152 

1090 

IH( I )«1 

15: 

CALL  ORDER!  MM1  ,  SIMP,  L,IH> 

152 

NITER=0 

1  53 

1100 

DONE=. TRUE. 

154 

NITER-HITEF+I 

155 

c 

1 56 

C 

COMPUTE  CENTROID  QF  THE  SIMPLE 

157 

C 

153 

CO  1110  1=1 ,MM1 

159 

1110 

CENTER! I  )=0.0 

1S0 

DO  1 130  1  =  1  ,MM 1 

16! 

IF:  I.EQ.  IH'.MMI  >  )  GO  TO  1130 

162 

00  1120  J  =  l  ,.MM 

163 

CENTER!  J  ) -CENTER!  J  >+5IllP<  J  ,1  ' 

MM+ROOT2  > 


A- 5 


3 


WS 

& 


% 

£ 

I 

V 

% 


I 


to 

ii 


$ 

u 

* 

$ 

£: 

» 

*3 


•JSEH 


164  1120  CONTINUE 

165  1130  CONTINUE 

166  C 

167  C  TRY  SPECULAR  REFLECTION. 

158  C 

1 69  QO  1140  1-1  ,MM I 

1 70  CENTER! I >«CENTERi I  )/MM 

171  RNEXT! I  )=( 1 .0+ALFA)*CENTER< I >-ALFA»SIMP( I  ,IH(MMI  )) 

172  1140  CONTINUE 

173  CALL  FUNC< RNEXT  .MMI  > 

174  C 

175  IF(RNEXT(MMI  '<  .LE  .SIMP!  NM1  ,L(HMI  )>)  SO  TO  1150 

176  60  TO  1170 

177  1150  CALL  NEUURfMMt , SIMP, RNEXT, IH) 

179  C 

179  C  TRY  AN  EXPANSION. 

180  C 

18!  00  i  160  1-1  .Mil 

182  RNEXT!  I  i  “GAMMA* SIMP!  I  )  :•  +  <  1  . 0-6AMMA  >«CENTER< 

1 83  1160  CONTINUE 

CALL  FUNCf RNEXT  , MMI > 

1S5  C 

186  I F  < RNEXT* MM  I  > .LE . SIMP! MMI  ,L(M111  >)) 

137  X  CALL  NEWURiMMl .SIMP .RNEXT , IH ) 

183  SO  TO  I2F0 

189  C 

190  1173  IF! RNEXT  (Mill  ) .  LE  •  SKIP!  AM  1  , IH( MMI  11)  60  TO  1190 


1SI 

GO  TO  11 90 

192 

1130 

CALL  NEUURi MMI .SIMP  .RNEXT . IH > 

193 

60  TO  12S0 

194 

C 

195 

C 

TRY  A  CONTRACTION. 

196 

c 

197 

1190 

00  1200  1  =  1  .MM 

198 

1230 

RNEXT < I )-BETA*SIMP(  I , IH<  MM1  ) )  +  < 1 .0-SETA  >«CENT£R<  I  ) 

199 

CALL  FUNC( RNEXT ,MM1 ) 

200 

C 

201 

IF «  RNEX  T(  Pim  S.LE.SIMPIMMI  ,IH(MM1  )!>  GOTO  1210 

202 

SO  TO  1220 

203 

1210 

CALL  NEU'.'R'.UM!  .SIMP  .RNEXT  ,IH  1 

204 

GO  TO  1250 

205 

2 

206 

Q 

DO  A  SHRINKAGE. 

207 

C 

209 

1220 

DO  1240  1*1  ,MM1 

209 

00  1230  J  =  1  ,MM 

2  10 

123C 

GIMP!  J  ,  I  3 1 IIP <  J  .1  :  »SIMP!  .1  ,L(NMl  ) )  >*SE1  A 

21  1 

1243 

CALL  FUNC (RNEXT. MM  1  ) 

212 

1250 

CONTINUE 

213 

1270 

FORMAT! 14/ > 

214 

I2S0 

FORMAT (612.4) 

215 

IFTLEQ)  LCNT =LCNT  + 1 

216 

IFC .MOT. LEO)  LCNT-0 

217 

OLOWL-OLOU 

218 

C 

219 

WRITE! LUQ  ,  1 270 )  NITER 

utj'ij'uru1 


223  1308 

224 

225 
225  C 
227  1305 


254 

255 

25S  132 


269  1420 

270 


00  1308  !•*)  ,MM1 

IFd.NE.KMI!  WRITE! LUO , 1 280 )  RNEXT( I  ) 

IF!  I.EQ.MM1  ,AND.RN£XT!MM1  J.LT.DLOW)  WRITE!  LUO  ,  I  280  )  RNEXT!  I  ) 
FINAL! I >»RNEXT<  I  ) 

WRITE!  LUO  ,1420  > 

IF< RMEXT ( MM1  ).LT.DL0W>  OLOW-RNEXT! Mill  > 

LEQ“! DLOUL . EQ.DLOU > 

CALL  OkOER!HM1  .SIMP.L.IH) 

00  1 320  J-t  ,MM 

ERRCF.i  J  >=(  SIMP'  J  ,IH(  I  )  >-  =  IHP 4  j  ,L(  J  )  )  >/SIMP!  J  ,IHi  J  >  ) 

CONTINUE 

WRITE!  LUO  ,♦  )  (ERROR!  J),J  =  1  ,MM  ) 

FOLLOWING  IS  PRESENT  AUTHORS’  CRITERION  FOR  CONVERGENCE , 

BASED  ON  EXPERIENCE  WITH  NONLINEAR  FUNCTIONS  IN  FITTING 
APPLICATIONS.  ESSENTIALLY,  THIS  CRITERION  SAYS  THAT  IF  NUMBER  OF 
CONSECUTIVE  UNSUCCESSFUL  ATTEMPTS  FIND  A  MEW  DIRECTION  AND 
SHAPE  TO  DECREASE  THE  VALUE  OF  THE  FUNCTION  REACHES  MM1»3, 

THEN  CALL  IT  A  GAME  AND  RETURN  TO  MAIN  PROGRAM.  THE  SIMPLEX  HAS 
EITHER  FOUND  A  MINIMUM  AND  IS  SPINNING  ITS  WHEELS,  OR  IT  IS  EATING 
UP  TIME  DUE  TO  POORLY  CHOSEN  INPUTS.  TRY  A6AIN  WITH  DIFFERENT 
INITIAL  ESTIMATES  AND  STEP  SIZES ,  IF  USER  SUSPECTS  SOMETHING 
UNUSUAL.  IN  FACT,  IT  IS  ALWAYS  GOOD  PROCEDURE  TO  VERIFY 
RESULTS  BY  RERUNNING  WITH  DIFFERENT  INPUTS,  AND  TO  INSPECT 
ITERATION  HISTORY  FILE  ITER. OUT. 

IF! ICON . ED .2  .AND.  LCNT.5E.MMt »3  !  GO  TO  1250 

THE  FOLLOWING  IS  A  DIFFERENT  TYPE  OF  CONVERGENCE  CRITERION. 

IF! ICON.EQ. 1  )  THEN 
K*0 

00  1322  J  =  1  ,MM 

IFi  ERROR! J  1 . LT . RMAXER' J  )  )  K-K  +  1 
CONTINUE 

IF(K.SE.MM)  GO  TO  1350 
END  IF 

IFtNITE® .6T.MAXITR  )  GO  TC  1350 
TRY  AGAIN. 

GO  TO  I  I  00 

CONVERGENCE  OBTAINED ,  OP  ITERATION  LIMIT  EXCEEDED,  ETC. 

CONTINUE 
WRITE! LUO, 1420) 

FORMAT! ’  ' ) 

RETURN 

END 


1 


s 


■ 

I 


SUBROUTINE  FUNC<  COEF  ,MM1 


272 

273  C 

274  C 

275  C 

276  C 

277  C 

278 

279 
288  C 
281 
282 
293 

284 

285 

286  C 

287  C 

288  C 

289  C 

290  C 

291  C 

292  C 

293  C 

294  999 


A  SIMPLE  TEST  FUNCTION,  WITH  THREE  UNKNOWN  COEFFICIENTS 
COEF(J),  J-1.3,  IS  TO  8E  MINIMIZED.  NOTE  THE  CURRENT  FUNCTION 
VALUE  (THE  MM+1  th  VERTEX  OF  THE  SIMPLEX)  IS  STORED  IN  COEF(MM+l  ) 


DIMENSION  COEF< | | ) 

IMPLICIT  DOUBLE  PRECISION  (A-H.Q-Z) 


TERMI  *=COEF(  !  '-10 
TERM2=CO£f (2  >-2.0 
T£RM3=COEF<3>-3.0 

COEFiMM! >-TERKt»*2  +  TERM2*»2  +  T£RM3»«2 
RETURN 


CONSTRAINT  CONDITION  EXCEEDED  -  SEND  BACK  A  MESSA6E  TO 
SUBROUTINE  SIMPLX.  BY  TAGGING  ON  A  LAR6E  VALUE  TO  THE  FUNCTION. 

NOTE  -  USER  MUST  PROTECT  AGAINST  SUCH  OCCUFRENCES  AS  TRYING 

TO  TAKE  LOGARITHMS  OF  0  OR  NEGATIVE  NUMBERS,  ETC.  TREAT  THEM 
AS  CONSTRAINTS,  AND  DIRECT  PROGRAM  TO  STATEMENT  999,  WHEN 
NECESSARY. 


CONTINUE 

COEFi MM  1  >*COEF( MM1  l+l .030 

RETURN 

END 


rraji  rjoywn 


39  SUBROUTINE  NEWVP.I KM1  ,S IMP  .RNEXT  . IH  ) 

30  C 

31  0  THIS  SUBROUTINE  STORES  THE  NEU  CHOICES  FOR  VERTICES  ,  WHICH  HAVE 

02  C  RESULTED  IN  A  DECREASE  IN  THE  VALUE  OF  THE  RESPONSE  FUNTION . 

03  C 

04  DIMENSION  SIMP! I !  ,1  >,RNEXT< I  >,IH(I > 

05  IMPLICIT  DOUBLE  PRECISION  IA-H.O-Z) 

06  0 

D7  DO  3013  1  =  1  ,MMI 

08  SIMP!  I  ,IH(MM:  )  )-RNEXT<  I  ) 

09  30! 3  CONTINUE 
1 0  RETURN 

1  1  ENO 


v  v.v  ,  v  v.v  ^vv.v  v 


a 


v3 

3 


312 

313 

314 

315 

316 

317 

318 

319 

320 
32! 

322  C 

323 

324 

325 

326 

327  4010 


330 


SUBROUTINE  ORDER! MM  .SIMP.L.IH? 


THIS  SUBROUTINE  DETERMINES  CURRENT  LOWEST  ANO  HIGHEST  VERTICES . 
A NO  STORES  ORDER  INFORMATION  U5IN6  SECOND  INCEX  OF  SIMP'J.I). 
THIS  WILL  BE  USED  3V  SUBROUTINE  SfMPLX  IN  DETERMINING  THE  NEXT 
TRIAL  SHAPE  AND  LOCATION  OF  THE  SIMPLEX. 


DIMENSION  SIMP!  I  1  ,1  )  ,Lt  I).IH(  1  1 
IMFUCIT  DOUBLE  PRECISION  (A-H.0-2) 


329  4020 
329 


DO  4020  J=> .MM! 

DO  4010  I-l .MMI 

IFiSIMPI  J,I  I.LT.SIMPI  J.U  J  )  J  )  LI  J  )-I 

IFISIMPl J  .1 >. ST, SIMP! J  ,IH( J  >>>  IH  < J  >“ I 

CONTINUE 

CONTINUE 

RETURN 

END 


3 


a 


V, 


1 


h 


0 


% 

>5 


Rl 


y. 


■:a 


> 


I 


k  n 


DISTRIBUTION 


IIT  Research  Institute 
ATTN:  GACIAC 
10  W.  35th  Street 
Chicago,  IL  60616 

US  Army  Materiel  System  Analysis  Activity 
ATTN:  AMXSY-MP 

Aberdeen  Proving  Ground,  MO  21005 

AMSMI-RD,  Dr.  McCorkle 
Dr.  Rhoades 

AMSMI-RD-CS-T,  Record  Copy 
AMSMI-RD-CS-R,  Reference 
AMSMI-RD-DE,  Jerry  D.  Smith 
AMSMI-GC-IP,  Mr.  Bush 


DIST- 1/ (DIST-2  Blank) 


